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BACKGROUND 


MGB Code New Modifications 



A number of modifications have been proposed in order to improve the jet noise 
prediction capabilities of the MGB code. This code which was developed at 
General Electric, employees the concept of acoustic analogy for the prediction of 
turbulent mixing noise. The source convection and also refraction of sound due to 
the shrouding effect of the mean flow are accounted for by incorporating the high 
frequency solution to Lilley's equation for cylindrical jets (Balsa and Mani). The 
broadband shock-associated noise is estimated using Harper-Bourne and Fisher's 
shock noise theory. The proposed modifications are aimed at improving the 
aerodynamic predictions (source/spectrum computations) and allowing for the non- 
axisymmetric effects in the jet plume and nozzle geometry (sound/flow interaction). 
In addition, recent advances in shock noise prediction as proposed by Tam can be 
employed to predict the shock-associated noise as an addition to the jet mixing 
noise when the flow is not perfectly expanded. Here we concentrate on the 
aerodynamic predictions using the PARC code with a k-e turbulence model and the 
ensuing turbulent mixing noise. The geometry under consideration is an 
axisymmetric convergent-divergent nozzle at its design operating conditions. 
Aerodynamic and acoustic computations are compared with data as well as 
predictions due to the original MGB model using Reichardt's aerodynamic theory. 
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MODELING APPROACH 


• Source Spectrum Calculations 

- Lighthill’s Acoustic Analogy 

- Ribner and Batchelor Assumptions 

- Calculation of Source Strength and its Spectrum Using CFD 

• Sound/Flow Interaction 

- High Frequency Asymptotic Solution to Lilley’s Equation for a Mul- 
tipole Source Convecting in an Axisymmetric Parallel Flow ( Balsa 
and Man!) 

- Calculation of Directivity Factors Based on CFD Results 
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GOVERNING EQUATIONS 


Lightbill’s Equation 


d ljL - <?v 2 p = JfUiL 

dt 2 P dxidij 


Tij = pViVj + 6ij(p - C 2 p) - e,j 
fdVi^dVj 2 dVk. 

6<i ~ ^dxj + dxi 3 6,1 dxj 

1 „ pUL 

Re-'— 

ldp- d p=(f s )ds 


• The effects of source convection and refraction are included in the source 
term. 
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LiUey’s Equation 


D ( D 2 a _ d jdo . dVj d 2 da 
Dt ( Dt 2 dxi * dxi ) + 2 dx i dxj ° dn 
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d d 
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The effects of source convection and refraction are included in the operator 
term of Lilley's eq. 
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Source Spectrum Calculations 


Mean-square sound pressure autocorrelation in the far field due to a finite 
volume of turbulence (in absence of convection and fluid shielding) 

d* 
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SOURCE SPECTRUM CALCULATIONS 


- Fourth-order velocity correlation tensor 

Sim = mWt = r mvi) dt 

J — oo 


- Source strength (Quasi- incompressible turbulence) 



- Reduction in order of correlation tensor (Ribner) 


Sij ki — Jb Sji -j~ S«S jk *+■ SijSkt 

/ +oo 

ViVjdt 

■OO 


- Separable second-order tensors 

•Sy(nl) = Rij (f)G(r) 


- Isotropic turbulence model of Batchelor 

&ij(d = re-^^*) 2 x{[ 1 - At/Lzffo+tMj/Ll} 

T=\wi, 

- Gaussian correlation time delay 

G(t) = e- (T/T - )S 

- Source spectrum component 

Ain(fi) ~ p 2 k*(Qr 0 ) 4 e ' •° > 

Jt 3 / 2 L 1 

Lz ~ T’ To ~ vl’ * = 2 ViVi 
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SOURCE SPECTRUM CALCULATIONS 


Characteristic time delay of correlation 

r ° ~ (du/dr) or To ~7 

Doppler shifted frequency 

Q = 27 r/C, M c = -5Af + PcMj 

C = \/(l - M c cos0) 2 + (a c fc' 5 /Coo) 2 


The proportionality constant in finding t 0 and the 
convection constants o: c and /3 C are dete rmin ed em- 
pirically. 


SOUND /FLOW INTERACTION 


Mean square pressure in the far-field 


j V 


+ 4a iy + 2o yy + 2oy Z )dy 


Source term 



~ (4irRCooCf(l - McosOf{ 1 - M c cos0) 2 

Shielding function 

2 (1 - Mcos^) 2 (% l ) 2 - cos 2 0 

** ^ ^ (1 - A/ C cos0) 2 

Af(r) = l/(r)/Coo, Af c = £/ c /Coo 

Directivity factors a XI ... are functions of the Shield- 
ing factor g 2 {r). 


23-7 




Geometry and Computational Grid 




Design Parmaters for C-D Nozzle 


Throat diameter 5.1 j n . 

Exit diameter 5.395 in. 

Distance from throat to exit 5.525 in. 

Exit velocity 2409 fps 

Ambient velocity 400 fps 

Pressure ratio 3.121 

Stagnation temperature 1 71 6 °R 
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MACH NUMBER 



Mach number contour plot for the convergent-divergent nozzle at the design 
condition (NPR = 3.12), using PARC code with a k-e turbulence model. The 
upstream conditions are specified at three diameters within the nozzle. The 
flowpath for this nozzle has been designed to obtain an isentropic, uniform and 
parallel flow at the exit for the design Mach number of 1 .4. 
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COMPARISON OF VELOCITY PROFILE 
WITH DATA 
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Comparison of velocity profiles with data at four different axial locations. The axial 
and radial dimensions are normalized with respect to the jet throat diameter D eq . 


23-10 





VELOCITY PROFILE 
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Comparison of velocity profiles with data on the (a) centerline (b) lip-line. 
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TURBULENT INTENSITY 



Turbulent intensity contour plot. The contour levels are normalized with respect to 
square of ambient sound speed. 
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COMPARISON OF PARC TURBULENT INTENSITY 

WITH DATA 
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Comparison of PARC turbulent intensity profiles with data and predictions due to 
Reichardt's theory on the (a) lip-line (b) X/D = 8.21. Radial distance Y is measured 
from the centerline and all percentages are based on the jet exit velocity Uj. Figure 
(a) shows that a maximum level of 13% is predicted by both prediction methods 
although the Reichardt's theory shows a much faster decay along the lip-line. The 
agreement between PARC and data is reasonably acceptable. 
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COMPARISON OF PARC TURBULENT INTENSITY 
WITH REICHARDT’S SOLUTION 



Comparison of turbulent intensity between (a) Reichardt's theory (b) PARC-ke. The 
radial profiles of turbulence can be compared at various axial locations. Figures (b) 
show that the centerline value of turbulent intensity peaks at X/D = 1 4.4 which is 
nearly twice the length of the potential core for this jet and decays farther along 
the jet. This is in agreement with experimental observations. Reichardt's model, 
on the other hand, fails to properly predict radial profiles of turbulence farther 
downstream of the potential core. 
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STROUHAL NUMBER 



The Strouhal number based on the inverse of the characteristic time delay of 
correlation nondimensionalized with respect to jet exit velocity and diameter. 
Figures (a) and (b) are concluded from PARC results using the two definitions el k 
and 3U/3r respectively. They show that, outside a proportionality constant, similar 
results can be obtained for the correlation time factor as a function of source 
location. Figure (c) is based on Reichardt's aerodynamic predictions and should be 
compared with Figure (b). 
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COMPARISON OF OVERALL SOUND PRESSURE 
LEVEL DIRECTIVITY WITH DATA 
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The overall sound pressure level directivity (OASPL) as estimated from PARC-ke is 
compared with data and predictions of Reichardt aerodynamic theory (on a 40 foot 
radius). The characteristic time delay r 0 is obtained from 1/r 0 = 2(e/k). The 
convection constants a c and 0 C are determined empirically. 
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COMPARISON OF SPECTRAL COMPONENTS 

OF NOISE WITH DATA 
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Comparison of noise spectra with data (based on 1 13 octave center frequency) at 
various observation angles. Band number 24 corresponds to 1 kHz. 
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COMPARISON OF SPECTRAL COMPONENTS 

OF NOISE WITH DATA 
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Comparison of noise spectra with data (based on 1/3 octave center frequency) at 
various observation angles. Band number 24 corresponds to 1 kHz. 
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SUMMARY 


Source Strength has been successfully pre- 
dicted using PARC code with a k-e turbulence 
model 

The limitation on aerodynamic grid selection 
has been removed by adopting a two-stage 
aerodynamic and acoustic algorithm 

The time-delay of correlation was calculated 
directly from kinetic energy of turbulence and 
its dissipation rate 

Constants used in supersonic convection factor 
need to be determined empirically 

The SPL directivity and spectra demonstrate 
good agreement with data especially at angles 
where mixing noise is dominant 
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